Reprinted from Journal of Applied Physics, Vol. 35, No. 6, 1764-1767, June 1964 
Copyright 1964 by the American Institute of Physics 




Printed in U. S. A. 


N 64 27760 


Superposition of Semi-Infinite Solenoids for Calculating 
Magnetic Fields of Thick Solenoids 
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A simple method is given for calculating the magnetic held components in or around any thick finite 
solenoid by superposition of fields of semi-infinite solenoids with zero inner radius. Equations and graphs 
are presented for the field components of such semi-infinite solenoids. From these graphs, the fields for 
solenoids of practical interest can be obtained with errors of less than a few percent. Greater accuracy is 
possible if numerical tables are used instead of graphs. 
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INTRODUCTION 

T HE magnetic field components off the axis of a 
thick solenoid cannot be calculated easily except 
by electronic computers. Tables of these fields, prepared 
from computer output, would often be useful, especially 
in coil design. Such tables would eliminate further use 
of computers in many cases. But the field components 
of a thick solenoid depend on four variables, which are 
the radial and axial coordinates of the field point and 
the two parameters a and /3 needed to specify the 
solenoid shape, where a — outer diameter/inner diameter 
and £= length/ inner diameter. Tabulating this function 
of four continuous variables would require an un- 
reasonable amount of space for complete coverage of 
the range of variables of common interest, although 
some tables and graphs exist which partially cover the 
ranges of the four variables. 1-4 

However, a comprehensive presentation in terms of 
only two variables is possible if a superposition method 
is used to find the desired field from tables or graphs of 
the field of a semi-infinite solenoid with zero inner 
radius. This semi-infinite solenoid is defined as an 
axially symmetric, uniform, azimuthal current density 
which extends from the axis of a cylindrical coordinate 
system out to R~b and from Z = 0 to Z= — oo (Fig. 1). 
It differs from ordinary solenoids in that it is semi- 
infinitely long and has no cylindrical hole. Each field 
component of this special solenoid can be expressed 
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Fig. 1. Semi-infinite solenoid with zero inner radius. Current 
density is uniform, extending from axis to R~ b and from Z— 0 
to Z — — » . Curved arrows indicate direction of current 
circulation. 
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nondimensionally, computed electronically, and tabu- 
lated or graphed as a function of only two variables, the 
nondimensional field -point coordinates, r—R/b and 
z—Z/b. Any desired finite solenoid is then treated as a 
superposition of four of these semi-infinite solenoids, 
and any field component is obtained by algebraic addi- 
tion of four numbers derived from the tables or graphs, 
each number corresponding to the contribution of one 
semi-infinite solenoid. 

SUPERPOSITION OF FOUR SEMI-INFINITE SOLE- 
NOIDS TO FORM A FINITE SOLENOID WITH 
FINITE THICKNESS 

Figure 2 shows how four semi-infinite solenoids of 
appropriate sizes and locations can be superimposed to 
form a thick, finite solenoid. Curved arrow r s show the 
directions of circulation of the currents. Cancellation of 
oppositely directed currents occurs in many regions, 
leaving only the desired finite solenoid. 

The field of the thick, finite solenoid can be expressed 
in terms of the fields of the four contributing semi- 
infinite solenoids in the following manner. Consider 



Fig. 2. Formation of finite solenoid from four semi-infinite 
solenoids. Curved arrows indicate directions of current circulation. 
Position of typical field point (X) is shown with respect to the 
desired finite solenoid and with respect to each semi -infinite 
solenoid. Note that since nondimensional tables or graphs are to 
be used, each semi -infinite solenoid effectively has its own coordi 
nate system, with its own radius as the unit of length. 
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either the axial or the radial component of the field of 
the infinite solenoid. To simplify notation, no symbols 
will be used to distinguish radial and axial field com- 
ponents; all equations in this section are valid for either 
component. Let H h H 2 , H 3 , and H 4 represent the values 
of the desired field component which would be con- 
tributed by the four semi-infinite solenoids marked 1, 
2, 3, and 4, respectively, in Fig. 2 if each had the con- 
ventional sense of current circulation. The actual fields 
contributed by the solenoids numbered 2 and 4 in the 
figure are thus —H 2 and ~H 4 , because the current 
densities have a reversed sense in these two semi- 
infinite solenoids. The field of the desired thick finite 
solenoid is simply the sum of the contributions of the 
four semi-infinite solenoids: 

H=Ih-H2+H z -H 4 . ( 1 ) 

If each term is divided by the current density J 
(assumed to be uniform) and by the inner radius s of 
the finite solenoid, a nondimensional expression results: 

H/Js = H \/ J s H 2 /Js+H,/Js~ H,/Js. 

Let the radii of the semi-infinite solenoids be called 
b h £ 2 , £ 3 , and £> 4 . Note that bi=b\=s and bi=b 2 =as. 
Then 

H!Js = aH l fJb l -aH 2 /Jb 2 +HJJb i -H,/Jbi. 

Let the nondimensional quantities such as Hi/Jbi be 
denoted by h h etc. Then 

H/Js=ahi—ah 2 - s rh3—h4. ( 2 ) 

The quantities k h h 2 , and in are the nondimensional 



Fig. 3. Nondimensional axial field of semi-infinite solenoid. 
Note: For z negative and r>l, no curves are shown; for such 
values, use the relation: h z (r , — z) = — h 2 (r,z ) for r> 1. 
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Fig. 4. Nondimen- 
sional radial field of 
semi-infinite solenoid. 
Note: Values for 2<0 
are not shown because 
Jlr (r, ~ z) =// r (r,z). 



fields of the four semi-infinite solenoids.. These non- 
dimensional fields can be approximately evaluated by 
using the graphs in Figs. 3 and 4 if the proper coordinate 
systems are used for locating the field point. The most 
convenient coordinate system to use with any particular 
semi-infinite solenoid is a cylindrical one with the origin 
at the center of the end face. It is convenient to non- 
dimensionalize the radial and axial coordinates by 
expressing them in tenns of the radius of the semi- 
infinite solenoid as a unit of measurement. Giving each 
semi-infinite solenoid its own origin of coordinates and 
unit of measurement means that the field point will have 
different coordinates with respect to each of the four 
semi-infinite solenoids. Figure 2 shows the coordinate 
system to be used for each semi-infinite solenoid : (ri,Zi) 
for solenoid 1, ( 72 , 22 ) for solenoid 2; etc. Note that the 
unit of length for measuring r\ and Z\ is b \ ; for r 3 and 23 
it is bz\ etc. Making these field point: coordinates 
explicit in Eq. (2) gives 

H/Js = ah (r i,z 1 ) — ah (r 2 ,z 2 ) + h (r- iy z 3) — h (r 4 ,z 4 ) . (3 ) 

Equation (3), containing only nondimensional terms 
which are expressed in terms of nondimensional co- 
ordinates, can be used to calculate the field components 
of any solenoid of finite length and thickness from 
h r (r,z) and h 2 (r y z). The functions h r (r,z) and h z (r,z) 
are the nondimensional radial and axial components of 
the magnetic field intensity of the semi -in finite solenoid, 
and their derivations are sketched at the end of the 
paper. They have been computed numerically and are 
plotted in Figs. 3 and 4. 

An example should clarify the procedure for calculat- 
ing a field by using Eq. (3) and the graphs. Suppose the 
axial field component is desired at the indicated field 
point for the finite solenoid shown in cross section in 
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EQUATORIAL PLANE} 



Fig. 5. Cross section of solenoid used in example calculation. 
Numbered corners correspond to locations of the edges of the end 
planes of the like-numbered semi-infinite solenoids of Fig. 2. 

Fig. 5. The inner radius of this finite solenoid is s. Each 
comer in the cross section is numbered and indicates the 
location of the circular edge of the end plane of the semi- 
infinite solenoid with the same number in Fig. 2. This 
numbering also corresponds to that of Eq. (3). Compar- 
ing Figs. 2 and 5, one can see that ri— 0.5/2 = 0.25, 
r 2 = 0.5/2= 0.25, r 3 = 0.5, r 4 = 0.5, Z\=— §=—0.5, 

2 2 =|= 1 . 5 j zz=3, and 2 4 = — 1. Since a= 2, Eq. (3) 
becomes 

H z /Js = 2h z (0. 25, - 0.5) - 2h z (0.25,1.5) 

+A,(0.5,3)-A,(0.5, -1). 

Figure 3 yields the necessary values of h zy hence 

HjJs = 2 X 0.62 - 2 X 0.030+ 0.008 - 0.45 = 0. 74. 

Although H/Js is a nondimensional quantity, the ex- 
pressions used in calculating h T (r,z ) and h„(r y z) were 
derived by using rationalized mks units, so a rationalized 
system of units must be used in these computations. 
Supposing /= 10 7 A/m 2 and 5= 0.05 m and noting that 
B=jjl 0 H in rationalized mks units results in the follow- 
ing for the axial component of magnetic induction : 

^= M t^X0.74-4TrX10- 7 X10 7 X0.05X0.74=0.46 

Wb/m 2 , or 4600 G. 

Reference 1 has a table of the fields of this particular 
solenoid and gives the more exact value of 0.7405 for 
H z /Js. The error in the above result is fortuitously 
much less than 1 %, but errors for other solenoids 
that are neither too short nor too thin will be typically 
only a few percent if the graphs presented herein are 
used. However, for very short solenoids or for very thin 
ones, two terms in Eq. (3) may be nearly equal but 
opposite in sign, causing a loss of one or more leading 
digits through subtraction. At the end of a calculation, 
it can easily be seen whether accuracy has been lost in 
this way. 

For calculating fields of such short or thin solenoids 
or to get greater accuracy for any solenoid, a tabulation 
of the field components of the semi-infinite solenoid is 
needed rather than graphs. Five-place tables and larger 
reproductions of Figs. 3 and 4 are to be published. How- 
ever, approximate field values can be obtained from 
Figs. 3 and 4 for a very wide range of the shape parame- 
ters a and j3. A remarkable feature of this method is 


that so much information about solenoids can be pre- 
sented on these graphs and can be handled without 
further use of a computer. 

The field of the semi-infinite solenoid has several 
properties which may be useful in calculations and some 
of which have been used to make Figs. 3 and 4 more 
compact. These are: 

(1) For z— 0 the nondimensional axial field com- 
ponent h z varies linearly from the value 0.5 on the axis 
to zero at r— 1. For z~ 0 and r> 1 the axial field is zero. 

(2) For a given value of r, the sum of the nondimen- 
sional axial field components at plus z and at minus z 
gives twice the value for 2=0 at the same r: 
kz(r,z)+h M (r, ~z) = 2h z (r y 0). 

(3) For r> 1 the nondimensional axial fields at posi- 
tive and negative values of 2 for a given radius are equal 
in magnitude but opposite in sign: h z (r,z)= — h z (r , —z) 
for r> 1. 

(4) For any given r the nondimensional radial field 
components are equal both in magnitude and in sign for 
positive and negative 2 : h r (r y z) — h r {r , — 2 ). 

CALCULATION OF THE FIELD COMPONENTS OF 
THE SEMI-INFINITE SOLENOID 

Let R , Z, and §> be the radial, axial, and azimuthal 
coordinates of the field point, as in Fig. 1, and let a, /, 
and 0 be the coordinates of the source point. The origin 
of the coordinate system lies at the center of the end 
plane of the semi-infinite solenoid. The azimuthal co- 
ordinates 0=0 and <F=0 are taken as coincident. The 
only nonvanishing component of the magnetic vector 
potential is the azimuthal component A$, which in 
rationalized mks units is 

(i<J f b f r a cosddd 

A^— / dll dal . 

2ir J — ao Jo Jo [_(Z-l) 2 A-R?A-a 2 —2aRcosdJ* 

Then by using H= (l//*o) curl A, expressions for H z and 
H t can be obtained. Let £=Z— l, then 

J r b r T a cosOdS 

H r =—\ da , 

2tJ 0 Jo [_Z 2 +R 2 +a 1 ~2aRcosey 

and 

J r 50 r b r * acosddS 

H z = — d£\ dal . 

2t tJ z Jo Jo i?+R 2 +a 2 -2aR cos0]* 

The expression for H r can be integrated analytically 
with respect to a, and the expression for H z can be inte- 
grated analytically with respect to both £ and a. If the 
resulting expressions are divided by Jb and if the non- 
dimensional field point variables r=R/b and z—Z/b are 
substituted into the expressions, the following results: 

H T 1 f* 

h r (r,z )= — = — / F cos Odd 
Jb 2ir J 0 

1 n r 1 — t cos0 q 

4 / r cos 2 0 sinh -1 \dd. 

2i rJ 0 L (s 2 + r 2 sin^) U 
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Hz 1 f r r r2(l-r cos#-HF)“| 

— = — —2 / In — \d0 

Jb 2A Jo L \z\ + (z>+r>)* J 
r * sin 2 $dd 

+zr 2 / 

Jo (l+r 2 -2rcos#)F 
rz f * r(l~r cos#) |z| q 

4 / sin# tan --1 \dd 

\z\Jo L rsin#F J 

+ !(i- f+ |i_ f | ) [ ) 


where F= (z 2 +r 2 +l~2r cos#)*. 

These expressions were integrated by a computer. 
The second integral in h r {r,z) is improper for 2 = 0 ; 
consequently, for 2=0 it was integrated by parts ana- 
lytically to avoid computer difficulties. 

Five-place tables (to be published) of h z {r,z) and 
h r (r,z) were automatically printed by the computer and 
Figs. 3 and 4 were plotted from the tables. 
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